Linear response theory of interacting topological insulators 
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Chiral surface states in topological insulators are robust against interactions, non-magnetic disor- 
der and localization, yet topology does not yield protection in transport. This work presents a theory 
of interacting topological insulators in an external electric field, starting from the quantum Liouville 
equation for the many-body density matrix. Out of equilibrium, topological insulators acquire a 
current-induced spin polarization. Electron-electron interactions renormalize the non-equilibrium 
spin polarization and charge conductivity, and disorder in turn enhances this renormalization by 
a factor of two. Topological insulator phenomenology remains intact in the presence of interac- 
tions out of equilibrium, and an exact correspondence exists between the mathematical frameworks 
necessary for the understanding of the interacting and non-interacting problems. 



I. INTRODUCTION 



The understanding of insulating behavior has been rev- 
olutionized by the landmark discovery of topological in- 
sulators (TI),I^1S1 "which are bulk band insulators with 
spin-orbit induced conducting states on the surface (3D) 
or edge (2D). These states are a manifestation of Z2 topo- 
logical order: topology guarantees the existence in equi- 
librium of a crossing of bands connecting time-reversal 
invariant momenta, which is robust against smooth time- 
reversal invariant perturbations such as non-magnetic 
disorder and electron-electron interactions. The surface 
states of 3D TI are described by a Rashba HamiltoniarP 
with Dirac-cone like dispersion, and are gapless and chi- 
ral, with a well-defined spin texture (spin-momentum 
locking.) They carry a n Berry phase, which protects 
against back-scattering and thus localization, and is asso- 
ciated with Klein tunneling, a half-quantized anomalous 
Hall effeclPand a giant Kerr effect .I^Nontrivial topology 
makes TI a platform for the observation of Majorana 
fermionaSl and for topological quantum computing.^ 

The rise of topological insulators is following a close 
parallel to the rise of graphene a short time ago. Three- 
dimensional topological insulators have grown from non- 
existence to a vastly developed mature field involving 
hundreds of researchers practically overnight. Within 
this time span, chiral surface states started out as a 
mere theoretical concept, were predicted to exist in sev- 
eral materials and were subsequently imaged.!^ Unlike 
graphene, the Hamiltonian of topological insulators is a 
function of the real spin, rather than a sublattice pseu- 
dospin degree of freedom. This implies that spin dynam- 
ics is qualitatively different from graphene. Moreover, 
the twofold valley degeneracy of graphene is not present 
in topological insulators. Despite the apparent similari- 
ties, the study of topological insulators is thus not a sim- 
ple matter of translating results known from graphene. 
Due to the dominant spin-orbit interaction, topological 
insulators are also qualitatively different from ordinary 
two-dimensional spin-orbit coupled semiconductors. 



The topological order present in TI is a result of one- 
particle physics. In light of this, we recall that electron- 
electron interactions modify the effective mass, heat ca- 
pacity, and ground state energy of solids, as well as 
the response of solids to external magnetic fields.liS! In 
fact, electron-electron interactions can lead to sponta- 
neous magnetism in itinerant electron systems. The best- 
known example of this effect is Pauli paramagnetism in 
interacting electron systems. It is known from Fermi 
liquid theory that the Pauli paramagnetic susceptibil- 
ity is enhanced by electron-electron interactions. This 
can be derived rigorously using various types of linear- 
response formalisms, such as diagrammatic Kubo linear 
response theory, the Keldysh kinetic equation formal- 
ism, or density-matrix formalisms based on the Liou- 
ville equation. Interaction effects in systems with s trong 
spin-orbit interactions have been studied in 2D TPDISl 
and 3D TI,^^]!!!! and previously in spin-orbit coupled 
semiconductors.'22n2S] j^ topological insulators the focus 
has been on phenomena in equilibrium and in the quan- 
tum Hall regime .33 

In the mean time, transport in topological insula- 
tors has made enormous strides recently.'^SJ In initial 
experimental efforts, it appeared impossible to identify 
any signatures whatsoever of the elusive surface states. 
Yet lately experimental work on transport in topolog- 
ical insulators has begun to advance at a brisk pace, 
and is without doubt entering its heyday, in the way 
ARPES and STM work did two years ago. A beauti- 
ful experimentpZI recently detected the topological surface 
states of Bi2Se3, in which Sb was partially substituted 
for Bi to reduce the bulk carrier density to lO^^cm"'^. At 
large magnetic fields the surface states were clearly seen, 
with Shubnikov-deHaas oscillations depending only on 
the perpendicular magnetic field, and oscillatory features 
growing with increasing field. Another work showed that 
carrier densities can be tuned over a wide range with a 
back gate.'SS A more recent experimental brcakthrouglP^ 
investigated surface transport in thin films of 612803 of 
thickness ~ lOnm, observing Landau levels that evolve 
continuously from electron-like to hole-like. In another 



breakthrough, Kim et al^^ studied Bi2Se3 surfaces in 
samples with thicknesses of < lOnm, using a gate elec- 
trode to remove bulk carriers entirely and take both sur- 
faces through the Dirac point simultaneously. Ambipolar 
transport was observed with with well-defined p and n re- 
gions, together with a minimum conductivity of the order 
of e^//i, reflecting the presence of electron and hole pud- 
dles. Exciting deve lopments in HgTe transport have also 
been reported.l^iEIl 

Due to spin-momentum locking, the charge current 
flowing on the surface of a TI is intimately linked to its 
spin polarization.'22! Firstly, it is evident that an out-of 
plane spin polarization can be generated by a magnetic 
field or magnetization. However, an in-plane magnetic 
field cannot generate an in-plane spin polarization for a 
Dirac cone: it simply shifts the origin of the cone and can 
be removed by a gauge transformation.'^Sl On the other 
hand, the combination of spin-momentum locking plus an 
electric field can be understood as a net effective mag- 
netic field, which is in the plane of the TI, and generates 
an in-plane spin polarization. 

This paper presents a study of the role of electron- 
electron interactions in topological insulators in an elec- 
tric field, and their effect on the spin polarization gener- 
ated electrically in the plane of the TI. A fundamen- 
tal question is whether basic TI phenomenology sur- 
vives interactions out of equilibrium. It is known that in 
transport topology only protects against 6acA;-scattcring. 
Topological protection stems from time reversal symme- 
try, whereas transport is inherently irreversible. There- 
fore robustness against electron-electron interactions in 
equilibrium does not translate into the same robustness 
in transport. I will demonstrate that the effect of in- 
teractions can be absorbed by a renormalization of the 
non-interacting charge conductivity and spin polariza- 
tion, and the response is qualitatively the same. Topo- 
logical insulator phenomenology therefore remains un- 
changed by electron-electron interactions in the steady 
state. 

A multiband matrix formulation is imperative to cap- 
ture interband dynamics and disorder effects, which give 
a nontrivial factor of two to the renormalization factor 
appearing in the charge current and spin polarization. 
This paper uses an alternative matrix formulation of lin- 
ear response theory, which contains the same physics as 
conventional approaches and is potentially more trans- 
parent, relying on the quantum Liouville equation in or- 
der to derive a kinetic equation for the density matrix. 
This theory was first discussed for graphene monolayerPSI 
and bilayers,"^ and was recently extended to topological 
insulators including the full scattering term to linear or- 
der in the impurity density.'^' Peculiarities of topological 
insulators, such as the absence of backscattering, which 
refiects the tt Berry phase and leads to Klein tunneling, 
are built into this theory in a transparent fashion. In this 
work, the formalism of Ref. 33, is extended to account for 
electron-electron interactions via a mean-field approach. 
Since transport in non-interacting systems was studied 



in that work, only minimal overlaps required for consis- 
tency have been retained in this article. It is assumed 
that T = so that electron-electron scattering is absent. 
The theory assumes spTp/h 3> 1, where ep is the Fermi 
energy, located in the surface conduction band, and r the 
momentum relaxation time. The ph ysics considered here 
is distinct from spin-Coulomb drag,'23EZl which requires 
electron-electron scattering, and from previous work on 
transport in non-interacting TI.I321 

Electron-electron interaction effects have also been 
studied in graphene transport .132 The interaction physics 
discussed here is to be distinguished from that of 
graphene, since, as stated above, graphene is a multival- 
ley system, and its Hamiltonian is a function of the pseu- 
dospin, due to the sublattice degree of freedom, rather 
than the real spin. It is also important to realize that the 
mean-field Hartree-Fock treatment of interactions is par- 
ticularly advantageous in topological insulators, because 
formulating a large- A^ renormalization group expansion is 
a challenging task. This is because, whereas in graphene 
the spin and valley degeneracies yield N — g^g^j —4, but 
in topological insulators N = gsgv — 1- 

The outline of this paper is as follows. In Sec. |IT] a 
general effective Hamiltonian for interacting systems is 
introduced. The dynamics of the density matrix in inter- 
acting systems are discussed in a mean-field formulation 
in Sec. |III[ Following that, the effective one-particle ki- 



netic equation is derived in Sec. IV This is then solved so 



as to obtain the correction to the conductivity and its en- 
hancement due to disorder, followed by a brief discussion 
of current TIs, a summary and conclusions. 



II. EFFECTIVE HAMILTONIAN FOR 

INTERACTING SYSTEMS 



The many-body Hamiltonian is 



cs 



al3 



a^-yS 



(1) 



The two-particle matrix element V^% g in a general basis 
{(/)„(r-)} is given by 



V^l,s = jd'r j d\'<i>l{r)<i>l{r')V'-{r~r')Ur)Ur'). 

(2) 
Hermiticity implies VajS-yS = ^*fSaB ^'^'^ identity of elec- 
trons VajS-yS — Vj3aS-f The antisymmetrized form is 



Val3j5 — ^ {Val3fS — VaPSj)- 



(3) 



I will consider henceforth the crystal momentum repre- 
sentation, where a = ks. The electron-electron interac- 
tion is taken to be explicitly of the Coulomb form. The 



many-body Hamiltonian is written as H^'^ + V^, where 

ks 



v^ 



(4) 



' ^ks 






The one-particle matrix element Hkk'ss' includes band 
structure terms and disorder. The matrix element v{q) — 
v{q) is given by (e,- is the relative permittivity) 



v{q) 



e 
2eoerq' 



(5) 



The real v{q) arises from Coulomb interaction matrix 
elements between plane waves. 



V, 



kiSi,k2S2^k^S3,k4S4 



— f5siS4 Ss2S3 f5fcl+fc2,fc3+fc4 ^(^2 " ^3 ) • 



(6) 
The term with fc2 — ky, is canceled by the positive back- 
ground of the lattice, so w(0) = 0. 

In TI in the random phase approximation (RPA), ab- 
breviating q = k' — k, one replaces viq) = v(q) — > -^^, 

£(?) 
where s{q) the dielectric function. The polarization func- 
tion is obtained by summing the lowest bubble diagram. 
At T = the long- wavelength limit of the dielectric func- 
tion is^^' 



eiq) = 1 



ATTeo^rA \ q 



(7) 



The polarization function was also calculated in Ref. [121 
The screened electron-electron Coulomb potential has 
the form 



v{q) e^ 1 

£(<?) 2eoer y/ k^ + k'^ - 2kk' cos jkk' 



2 



(8) 



The Wigner-Seitz radius Vg (alternatively, the effective 
fine structure constant), which parametrizes the relative 



strength of the kinetic energy and electron-electron inter- 
actions, is a constant for the Rashba-Dirac Hamiltonian, 
given by r^ — e^ H^-KtQt^A). In addition to the electron- 
electron Coulomb potential, the matrix element IJ^k' of a 
screened Coulomb potential between plane waves, which 
will be relevant in transport below, is given by 



Ukk' — 



1 



Ze' 



2eoer |fe — fc'l + kTF 



(9) 



where Z is the ionic charge (which I will assume for sim- 
plicity to be Z = 1) and k^^p — kpr g/2 is the Thomas- 
Fermi wave vector, with kp the Fermi wave vector. 



III. DENSITY MATRIX 

The many-particle density matrix F obeys^ 

f + ilff.Fl^O^ (10) 

The one-particle reduced density matrix p is the trace 



POi = Tr(4cc^) = (^^c) = {P)-^e 



The reduced density matrix satisfies 



(11) 



^^C.+^[^le^^j^^^_^^[^.e^^t,^j^^O_ (12) 



dt 



In terms of the antisymmetric Coulomb matrix element 
Vap^s defined above, the last term on the LHS 

Q/37 



The many-electron average is evaluated as follows 



(13) 



t t 



a/3'y 
(c'lc'fjC^Cs) = {clcs){c'pCj) - {clc^){c'i^Cs) +Ga|3■y5■ 



{Ue.) 



(14b) 



The focus of this work is on the first two terms on the RHS of Eq. ( 14b ) , which represent the Hartree-Fock mean-field 



part of the electron- electron interaction. The remainder, Gap-yS, gives the electron-electron scattering term in the 
kinetic equation,™ is second-order in the interaction and vanishes at T = 0. I will treat the case of zero temperature 
and reserve electron-electron scattering for a forthcoming publication. To evaluate the Hartree-Fock factorization, 
note that Vap-ynP-fP cancels, and the remainder becomes 



{[Vee, cljC,;]) = 2^ [Va/Sjri (PCaP-yP ~ P'laP^fi) + Vfi^aC {Pir,PSia " PSir]P-,a)\ ■ 



(15) 



I will introduce two mean field terms by letting 
14/3777^7,3 = Aaif and 14/3777 P7Q = S^lf , then 

([Ke,ctcd) = [p,^^-^^-S*^nC.- (16) 



The effective kinetic equation becomes 



dp 
di 



[Hie,p] + j^[A 



MF 



6*^^,p] = 0. 



(17) 



The one-particle Hamiltonian is renormalized by 

Hllf = A''^ - B''^. (18) 

I emphasize that in the final analysis one is inter- 
ested only in the impurity average of p in the crystal- 
momentum representation. In general one may al- 
ways write pkk' = fkSkk' + gkJ^, where the fc-off- 
diagonal part, gkk', is eventually integrated out to yield 
the scattering term in any desired approximation. In 
the impurity average of Eq. (21), out of the commu- 
tator {[Himp),p])^,j only the terms [Himf)J] and 
[H^lf{g),g] survive, where to first order in the electric 
field Hgf{p)p = H!ff{po)pE + mlJ{pE)po- This im- 
plies that Hllf{g)g = Hlif{go)gE + Hlif{gE)go- In 
linear response 50 = 0, and we are left with H^l-^{f)f. 
Specializing to this term, spin indices are omitted and fk 
is treated henceforth as a 2 x 2 matrix in spin space. 

To determine H^l^ [f), we evaluate the two mean field 
terms. Beginning with A , with summation implied 
over repeated indices, 



Note that B^'^^ can be interpreted as an effective mag- 
netic field due to the Hartree-Fock mean field electron- 
electron interaction. This result reproduces the cor- 
rect exchange energy,GSj and yields exchange enhancement 
of Zeeman field-induced spin polarizations, as found in 
Fermi liquid theory. It is similar in spirit to the treat- 
ment of Ref. liTl 



Equation (12 1 is reduced to 



dp 
di 



H^'.p] 



\B 



MF 



P]- 



(21) 






VaP^nf-iP = Vk^s^M 



,k4^S4^Jk2 



— V{k2 — fca) ^SiS4 ^S2S3 Ski+k2M+k4 Ik^ 

.-. -> 0, since v{Q) = 0. 



(19) 



Therefore A^^^ vanishes in the most general case. Next, 
qMF jg given by 



^frf^ = VaM-nPia =^v{k - k' ) fk' 



(20) 



The single-particle Hamiltonian is renormalized by 
B^^^{f), which is itself a function of the single-particle 
density matrix. 



At this stage one may include explicitly disorder and 
driving electric fields in the one-particle Hamiltonian and 
write Hll, = H^kSkk' + HEkk' + Ukk' , where iJofe is 
the band Hamiltonian, HEkk' the electrostatic potential 
due to the driving electric field and Ukk' the disorder 
potential. The effective single-particle kinetic equation 
takes the form 



^ + ^ [H,k,fk\ + J{fk) = -^ [HE,,fk\ + ^ [Bi'^Jk]. 

I 



(22) 



One writes fk = fok + lEk + ft" . where /ofc is the equilib- 
rium density matrix, fEk is induced by the electric field, 
and ft"^ by electron-electron interactiona^Sl. 



Equation (22) is solved iteratively in B^'^ 

, fk 



Let the 
The approach is 



bare driving term Dk — ~i [HEk 
to solve the kinetic equation first with Dk as the source 
term. This will give a spin polarization. The spin polar- 
ization will give a nonzero Bk , which in turn will give 
an additional source term, referred to as Vf!^ in the next 
section. Then one solves the kinetic equation again with 



Vjf as the source term 

—T, 1" J [Hok, fEk] + J (/fife) = —-r [Hsk, fok] 



dt 



(23) 



[i?ofc,/r] + ^(/r) = i[^f^Jofe] 



h 



On the RHS of the second equation only the equilibrium 
density matrix /ofc appears because B^^^ is first order in 
the electric field. The iteration is continued to all orders 
in the Wigner-Seitz radius r^ (that is, to all orders in the 
effective fine structure constant.) 

We recall that electron-electron and electron-impurity 
potentials are screened, with screening treated in the 
random-phase approximation. The density-matrix for- 
malism used here is thus equivalent to the GW approx- 



imation. In the non-equilibrium diagram technique, the 
correction discussed in this work is obtained by includ- 
ing the real part of the Green's function due to electron- 
electron interactions.HS 



IV. KINETIC EQUATION FOR INTERACTING 

TI 

Henceforth I specialize to TI. The band Hamiltonian 
h 2Ak - 
Hok = -^ "■•^fc, where ftk — 7— S, with 9 the tangen- 
tial unit vector in polar coordinates in reciprocal space. 
Interaction with the electric field is given by He kk' = 

d 
(eE ■ f)fcfc'l = ieE ■ —- d(k - k') 1, with 1 the identity 

ak 
matrix in spin space. Uncorrelated impurities located at 

Ri are represented by Ukk' = Ukk'y^'^'^''''''^^' ' '^^^^ 

Ukk' the Fourier transform of the potential of a single 
impurity. I will write fk = n^l -t- S'fc, with Uk the num- 



The solution to Eq. (28) requires certain approxima- 



ber density and Sk = -^ Sk ■ cr the spin density. One 
decomposes Sk — Sk\\ + Sk_L^ where [HQk,Sk\\\ = and 
Sk\\ is the fraction of carriers in eigenstates of -ffofcj while 
Ski. represents interband dynamics, i.e. Zitterbewegung. 
Further, Sk\\ = {l/2)sk\\<Tk\\ and Ski_ = {1/2) Sk±(Jk±, 
with the matrices ak\\ — —c ■ 6 and ak± — cr ■ k. 



A. Single-particle kinetic equation 



The general single-particle kinetic equation is 



dSk\\ 
dt 



P\\J{Sk 



V. 



(24a) 



7 Q 

-^ + \ [Hk, Ski_] + Pj{Sk) = Vk^, (24b) 



where the driving term Vk 



eE dpi 



Ok 



and pok is 



the equilibrium density matrix. This equation is solved 
as an expansion in the small parameter h/{eFT), where 
the momentum relaxation time t is defined below. The 
leading-order term in this expansion is ex [h/{epT)Y^^' 
and is found from 



tions. With respect to the scattering potential one ex- 
pands in the small parameter h/{eFT). In the steady 
state in the Born approximation the leading term in the 
solution to the kinetic equation is oc [h/{eFT)]'--~^\ It is 
trivial to check that at finite doping the next term in the 
expansion, i.e. ex [h/{epT)Y^' , vanishes identically, which 
was demonstrated in Ref [551 A term ex [/i/(ei?T)]'") 
would appear in the weak localization regime, yet this 
correction is not relevant in the regime spT/h ^ 1 con- 
sidered in this work. 

The Born-approximation scattering term has the form 



Jih) 



h^\ 



dt'[U,e'^[Uj]e' 



kk 



(26) 

with J = 6' — 6 the angle between the incident and scat- 
tered wave vectors, k and k' respectively, and (()) denot- 
ing the average over impurity configurations. The projec- 
tions of J{fk) needed in this work have been determined 
befor^ 



P\\JiSk\\) = g^^^" / d'y\Ukk'\'^{sk\\ -Sfc'ji)(l + cos7) 

^_L-^('S'fc||) = '^^'^ j d-i\Ukk'?{sk\\ -.Sfc'||)sin7 

krii (Jk\\ /" I - i2 /- \ 

P\\J[Ski_) ^ c^t A id-i\Ukk'\ (sfe_L + Sfe'_L) sin7, 

(27) 
where 7 = 6*' — 6' is the angle between the incident 
and scattered wave vectors. The small Sk\\ and Sk_L are 
scalars, s^n = —Sk ■ and Ski_ = Sk ■ k. 

The scattering terms contain factors of (1 + C0S7) (re- 
fiecting the tt Berry phase) or sin 7, both of which pro- 
hibit backscattering and give rise to Klein tunneling. 
Since the current operator ex cr, only Sk is needed. In the 
absence of scalar terms in the Hamiltonian, J{fk) does 
not couple rik with Sk , and Eq. ( 22 ) makes evident the 



P||J(5fc||)=2?fc|| 



(25) 



fact that the interaction term does not couple Uk and Sk, 
thus Uk may be dispensed with for the remainder of this 
work. The equation satisfied by Sk is 



dSk 
dt 



^ [Hok,Sk] 



J{S^)^-'-[HE„Sk] + j^[B^ 



MF,{1) 
k 



(28) 



r 



With respect to the electron-electron interaction one also 
needs to define a perturbation expansion in order to solve 
Eq. (28), which is done in what follows. Within the ap- 



proximations used in this paper, this expansion can be 
summed exactly. The method of solution is summarized 
as follows. The kinetic equation first with B^'^ set to 



zero. This solution is already knowrPSI and gives a spin 
polarization, which in turn generates a nonzero B^^^, 
which itself yields a new driving term, and so forth. The 
full solution is found as a perturbation expansion in the 
electron-electron interaction, which can be summed ex- 
actly. 

To obtain the solution in the interacting case, it 
is therefore first necessary to solve the non-interacting 
problem. In the absence of interactionP^l the steady- 
state solution to the density matrix in the Born approx- 
imation is'^'^ 



S 



TeE-k a/o4 



Ek\\ 
1 

r 



Ah 

krii 
4J1A 



dk 



■^k\\ 



2^ 



(29) 



\Ukk'\^sm'^l- 



Above Tii is the impurity density, while the factor of sin 7 
represents the product (1 -t- cos7)(l — cos 7). The first 
term in this product is characteristic of TI and ensures 
backscattering is suppressed, while the second term is 
characteristic of transport, eliminating the effect of small- 
angle scattering. In non-interacting TI the Zitterbewe- 
gung contribution to the conductivity/spin-density (i.e. 

due to 5'^'fcj^ ) vanishes identically at finite doping. But 
in the interacting case it is necessary to consider both the 
electron and the hole bands to capture the spin dynamics. 
The solution obtained, Ssk = SEk\\: is fed into B^^^, 
which in turn generates a new driving term in the ki- 
netic equation. Each term in this expansion by the index 



a, thus Bi, 



The solution found in Eq. 



29 



sponds to a = 0, that is, in the non-interactmg case 
S 



Ek\\ 



,(0) 



'Ek\\ 



. The driving term due to electron-electron 
interactions is generically denoted V'^'^"'' . The decom- 
position Vj^j^"' — (1/2) d'^^°'' a ki_ is also used. The so- 
lution to the spin part of the density matrix to order a is 



,(") 



(") 



denoted by 5*^^. . The driving term 2?^, is always or- 
thogonal to Hok, therefore V'i^'^°'' 



V 



fie, {a} 
k± 



The kinetic 



equation for the solution S^'^ in the presence of electron- 
electron interactions can be written for each order as 



dS\ 



;e,(a) 
k± 



dt 



^ [Hk,Slj^ 



^ V 



ee,{a) 
k± 



(30a) 



PwJist 



e,(a)N 



-Pii j(5r/"^ (30b) 



where in Eq. (30b I it is understood that the RHS, found 



from Eq. (30a), acts as the source for the LHS. The scat- 



tering term does not appear in Eq. (30a) since, as was 
argued above, T^'j^'u = always. 



I will dwell first on the solution of Eq. ( 22 ) due to 



iJi, , i.e. first order in the interaction, which requires 



^k 

~,ee 

'-Efcll 



S^^i^^^ FromEq. (|29|), 



B 



MFXi) _ e^fci? 



ATTener 



1 
dl'l' 



^d-y 



S' 



e,(0) 
Ek' 



oCr Jo Jo 27r y/p^Y 



211' cos 7 -I- ^ 
(31) 



where / = (k/kp)- The term in B^ in which Sk 

Sok' gives a vanishing contribution. For E \\ x 



B 



MF,{1) 
k 



I^l\l,r,) = 



l^TrheQe, 
1 f d-f 



(cos 27 — 1) 



2vr [y/l + P -21 C0S7 + 



(32) 

In Ic {l,rg) the sign of the cosine term is flipped. Al- 
though Bj^ itself has a part ex cr^.n, this part drops 
out of the driving term in Eq. ( 33 ) , because one is work- 



ing to first order in the electric field and the commu- 
tator [B^'^'^^\Sk] -^ [B^^ ,Sok], and Sok « crk\\- 
The effective electron-electron interaction Hamiltonian 
B/^ therefore contributes a driving term orthogonal 

to i?ofej yielding a correction S'^^\^ to the density ma- 
trix. The scattering term does not appear in the equation 



for St, ■ Equation (|28|) takes the simple form 



dS 



.(1) 

Ek± 



rf, ^ f^iHo,S-,^l^]^'-[Br^^'\Sok]. (33) 
This equation is solved using the time evolution operator 



(1)/ 



,(1) _ eE^rsTfee {l,rs) 



'Ek± 



\mk 



/o+ singer • k. 



Another contribution stems from the projection 



P\\JK 



ee,(l). 



-P\\J\S\ 



fe_L 



,(l)l 



(34) 



(35) 



It is understood that the RHS, found from Eq. (34), acts 
as the source for the LHS. Straightforwardly 



Si 



,(1) 



d^( 
\mk 



eE-j^r^Tl^e {l,rs) 



/o+ cos 6 a ■ 9. 



(36) 



•^BfcV ^^"^ '^fcii contribute equally to the charge cur- 
rent determined below. In effect, scattering from S'^^j^j^ 
into S'^y doubles the contribution to the electrical con- 
ductivity due to S*^^^^ . 

eA 



The longitudinal current density operator j 

3 equivalent to a spin pc 
J. of the non-interacting 

- jl33. The first-order conductivity 
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the current density is equivalent to a spin polarization. 
The conductivity cr^^ of the non-interacting system is 
V\ fAkpT^ 

correction in the electron-electron interaction is 
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(37) 



where lie (j's) = / dl I)^J [l , r s) . The electrical current 

. . ."'0 . 
and non-equilibrium spin polarization are renormalized 

by electron-electron interactions. 



Equation (371 has been obtained to first order in the 
(screened) interaction. The source term due to (f^j} 
contains only e****, identical in structure to the non- 
interacting problenpS One solves for all higher terms 



in H^'^ by iterating steps (31)-(36), obtaining the exact 



result for the conductivity (and spin polarization) 



The general formula for the dimensionless integral /, 
for n > 1 is 
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In 2D v{q) oc 1/q, 
independent, and the 



while in TI r^ is density- 
screened Coulomb potential 



v{q)/e{q) oc l/q also. Thus H^'' does not introduce den- 
sity dependence: at larger densities the Coulomb inter- 
action is weaker. Solving for S'^p.j^ introduces a factor of 
l/rjfc, which is canceled by k in the 2D volume element. 
Thus, 2D physics and TI linear dispersion combine to 
ensure the renormalization is density independent. 

The renormalization reflects the interplay of spin- 
momentum locking and many-body correlations. A spin 
at k feels the effect of two competing interactions. The 
Coulomb interaction between Bloch electrons with k and 
fc' tends to align a spin at k with the spin at fc', equiv- 
alent to a ^-rotation - hence the driving term in Eq. 
( 33 1 is ex CTz . The total mean-field interaction tends to 



align the spin at k with the sum of all spins at all fc', 
and H^*^ encapsulates the amount by which the spin at 
fc is tilted as a result of the mean-field interaction with 
all other spins on the Fermi surface. The effective field 
rife tends to align the spin with itself. As a result of 
this latter fact, out of equilibrium, an electrically-induced 
spin polarization is already found in the non-interacting 
systerrPH. Given that the spins at fc and fc' are in the 
plane, interactions tilt the spin at fe in the direction of 
the spin at fc'. Thus far the argument helps one under- 
stand why, if there is no spin polarization to start with, 
electron-electron interactions do not give rise to a spin 
polarization. The mean-field result is zero, so there is 
no overall tilt on any one spin due to the spins of the 
remaining electrons. Interactions tend to align electron 
spins in the direction of the existing polarization. The 
effective ^-rotation explains the counterintuitive obser- 
vation that the renormalization is related to interband 
dynamics, originating as it does in Sk±- Many-body in- 
teractions give an effective fc-dependent magnetic field 
II i, such that for .E || x the spins Sy and —Sy are rotated 
in opposite directions. Due to spin-momentum locking, 
a tilt in the spin becomes a tilt in the wave vector: spin 
dynamics create a feedback effect on charge transport, 
renormalizing the conductivity. This feedback effect is 




0.2 0.4 0.6 0.8 

Wigner-Seitz radius rs 

Figure 1: Fractional change in the conductivity \5a^%/a^^\ for 
< Ts < 1. The current generation of topological insulators 
has Ts <C 1, so the theory presented in this work is applicable 
to these materials. 



even clearer in the fact that the projection —P\\J{Sk±) 
doubles the renormalization. This doubling is valid for 
any elastic scattering. 

I will discuss next the magnitude of this renormal- 
ization in currently known topological insulators. Sev- 
eral materials were predicted to be topological insu- 
lators in three dimensions. The first was the alloy 
Bii^a^Sba: !'^^*'*'^* followed by the tetradymite semiconduc- 
tors Bi2Se3, Bi2Te3 and Sb2 Tea .'^ These materials have 
a rhombohedral structure composed of quintuple layers 
oriented perpendicular to the trigonal c-axis. The cova- 
lent bonding within each quintuple layer is much stronger 
than weak van der Waals forces bonding neighboring lay- 
ers. The semiconducting gap is approximately 0.3 eV, 
and the TI states are present along the (111) direction. In 
particular Bi2Se3 and Bi2Te3 have long been known from 
thermoelectric transport as displaying sizable Peltier and 
Seebeck effects, and their high quality has ensured their 



place at the forefront of experimental attention^ Initial 
predictions of the existence of chiral surface states were 
confirmed by first principles studies of Bi2Se3, Bi2Te3, 
and Sb2Te3.HSl In the current generation of topological 
insulators, r^ is small. Currently e^ ranges between 30 
and 100 (200 for Bi2Te3), making r^ between 0.14 and 
0.46. The theoretical treatment adopted in this work 
is therefore applicable, and interactions provide a correc- 
tion to the steady-state response. A plot of \Sa^%/a'^^ \ for 
< Ts < 1 is shown in Fig. [l] from which it emerges that 
interactions may account for up to ~ 15% of the observed 
conductivity of surface states in the regime studied. I 
note that Heusler alloys were recently predicted to have 
topological surface states,S£las well as chalcopyrites,S3yet 
more work is needed to establish the size of r^ in these 
materials. 

At this stage in topological insulator research, the re- 
sults found in this work are interesting for conceptual 
reasons, since they demonstrate that TI phenomenology 
is unchanged by interactions. The electrical conductiv- 
ity/spin polarization has the same form as in the non- 
interacting case, with a renormalization that can be in- 
corporated into a redefinition of the spin-orbit constant, 
or alternatively of the Fermi velocity, and thus the den- 
sity of states. For large r^ a non-perturbative treatment 
that goes beyond the random phase approximation is nec- 
essary, yet such a theory must await materials progress. 
In this context, I would like to note that the growth of 
new materials is a nontrivial issue, and obtaining high- 
quality samples where only the surface electrons can be 
accessed in transport has proved to be a difficult task. 
It is especially important to recall that future work may 
initially be hampered by factors such as roughness and 
dirt inherent in solid-state interfaces. In addition, it re- 



mains true that all current TI materials are effectively 
bulk metals because of their large unintentional doping 
- at present, bulk carriers are only removed temporar- 
ily by gating. Discussing TI surface transport in such 
bulk-doped TI materials retains some ambiguity, since 
it necessarily involves complex data fitting and a series 
of assumptions required by the necessity of distinguish- 
ing bulk versus surface transport contributions. Real 
progress is expected when surface TI transport can be 
carried out unambiguously, without any complications 
arising from the (more dominant) bulk transport chan- 
nel. The immediate tasks facing experimentalists are get- 
ting the chemical potential in the gap without the aid of 
a gate, further experimental studies confirming ambipo- 
lar transport, and the measurement of a spin-polarized 
current. 



V. CONCLUSIONS 

I have demonstrated that, from the point of view of the 
non-equilibrium spin polarizations and charge current, 
TI behavior remains intact in the presence of interac- 
tions with only quantitative modifications. The conduc- 
tivity and spin polarization are renormalized by electron- 
electron interactions entering through a combination of 
intcrband dynamics and scattering. 
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